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We update and extend our previous work reconstructing the potential of a quintessence field from 
current observational data. We extend the cosmological dataset to include new supernova data, 
plus information from the cosmic microwave background and from baryon acoustic oscillations. We 
extend the modelling by considering Pade approximant expansions as well as Taylor series, and by 
using observations to assess the viability of the tracker hypothesis. We find that parameter con- 
straints have improved by a factor of two, with a strengthening of the preference of the cosmological 
constant over evolving quintessence models. Present data show some signs, though inconclusive, of 
favouring tracker models over non-tracker models under our assumptions. 



PACS numbers: 98.80.-k 



I. INTRODUCTION 



The nature of dark energy in our Universe remains 
unknown, and is likely to be the subject of intense ob- 
servational attention over the coming decade While 
a pure cosmological constant remains the simplest in- 
terpretation of present data, a leading alternative pos- 
sibility is the quintessence paradigm, whereby the ob- 
served acceleration is driven by the potential energy of 
a single canonically-normalized scalar field Q (for ex- 
tensive reviews of dark energy see Ref. Q). In this pa- 
per, we work under the assumption that quintessence is 
a valid description of observational data (an assumption 
to be tested separately) , and seek to impose optimal con- 
straints on the model via exact numerical computation. 
Our work provides an implementation of quintessence 
potential reconstruction, a subject developed by several 
authors [J, [f| [y], and by assuming a particular physi- 
cal model for dark energy is distinct from parameterized 
equation of state methods for reconstructing dark energy. 

In a previous paper we carried out a direct re- 
construction of the quintessence potential based on the 
supernova type la (SNIa) luminosity-redshift measure- 
ments made/collated by Riess et al. Q. The present 
paper updates and extends that work in three ways: 

1. We include additional data coming from cosmic mi- 
crowave background (CMB) anisotropics [8[ and 
baryon acoustic oscillations 0, as well as using 
newer supernova data from the SuperNova Legacy 
Survey (SNLS) [l(|. We do not use constraints 
from the growth rate of structure, which are not 
yet competitive with the data we do use. 

2. Where previously we approximated the quint- 
essence potential via a Taylor series, we now ad- 
ditionally explore use of Pade approximant expan- 
sions in order to test robustness under choice of 
expansion. 
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3. By studying the dynamical properties of models 
permitted by the data, we assess whether current 
observations favour or disfavour the hypothesis that 
the quintessence field is of tracker form, hence po- 
tentially addressing the coincidence problem. 

As we were completing this paper, a closely-related pa- 
per was submitted by Huterer and Peiris [ll|, who also 
reconstruct quintessence potentials from a similar compi- 
lation of current data. Although phrased in the language 
of flow equations, their approach, like ours here and in 
Ref. [f|, amounts to fitting the coefficients of a Taylor 
expansion of the potential. They do not consider Pade 
approximants. Their approach implies different priors 
for the parameters than the ones used in this paper, and 
they treat the scalar field velocity a little differently. Our 
results appear in good agreement, in particular our deter- 
mination that present data mildly favour tracker models 
over non-tracker models concurring with their conclusion 
that freezing models are mildly preferred to thawing ones 
(in the terminology of Ref. [l2|). 



II. FORMALISM 

A. Cosmological model 

We quickly review the set-up of Ref. [6( , which is con- 
ceptually straightforward. We assume that the quint- 
essence field <j) has a potential V (4>), which we expand in 
a series about the present value of the field that is taken 
(without loss of generality) to be zero. The quintessence 
field obeys the equation 



dV 



(1) 



with the Hubble parameter H given by the Friedmann 
equation 



8nG 



(Pm + Pt) ■ 



(2) 



2 



Here p m is the matter density and p^ — (j> 2 /2 + V{4>) 
the quintessence density. We assume spatial flatness 
throughout (as motivated by CMB measurements and 
the inflationary paradigm), though the generalization to 
the non-flat case would be straightforward. Since then 
f2 m + Qtj, = 1 we have the present boundary condition 



ifo = ±72 [(1 - O m )p c , - V(<h)] , 



(3) 



where subscript '0' indicates present value, and p c is 
the critical density. An important quantity, which de- 
termines the cosmological effects we consider from the 
quintessence field, is the equation of state 



V4> _ 4> 2 /2 - V{<j>) 



P4> 072 + V{4>) 
The priors we assume for our cosmology are 

"total = 1 i 

> 0, 

^kin < 1 , 

f^kinO > 1) < 0.5. 



(4) 



(5) 
(6) 
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where fikin = 8ftG(j. > /6H is the fraction of critical en- 
ergy density in field kinetic energy density. The last con- 
dition is a means of encoding that the field should not 
interfere too much with structure formation (as we do not 
use data sensitive to that) , and is discussed further in our 
previous paper @. The constraint on fikin is in practice 
applied up to the highest redshift for which we have data 
points, i.e. using CMB information z = 1089. When we 
use supernova data only, the upper limit is z — 2, as in 
our previous study. 



B. Parameterizations and priors 

To explore the space of potentials, we need to assume 
some functional form for the potential. We choose two 
classes of expansions, a Taylor series, and a Pade series, 
to parameterize the potential function V ((/>). In the ab- 
sence of a theoretical bias for the functional form of the 
potential, these expansions seem suitably general and 
simple to provide a reasonably fair sampling of the space 
of potential functions. 



1. Taylor series 

As in our previous study, we use a Taylor series to 
model the potential V(<p) as 



V(<t>) ^V + V 1 4> + V 2 4> 2 + 



(9) 



where <p is in units of the reduced Planck mass Mp with 
<j) presently zero. We will refer to a constant potential 



with non-zero kinetic energy allowed as a skater model, 
after Linder in Ref. [l3j . 

We put the following flat priors on the parameters: 



Vo>0, I Vx| < 2 , |V 2 |<5. 



(10) 



These priors are irrelevant for parameter estimation, as 
they are significantly broader than the high-likelihood 
region (this also applies to the corresponding priors for 
Pade series below). However, to assess how favoured 
tracker behaviour is, we do need to put some limits, so 
that we can sample a finite region of the prior parameter 
space (see further in Section [TV C 2p . 



2. Pade series 

In addition to the Taylor series expansion, in this paper 
we also use Pade approximant expansions in order to 
test the robustness of results to the method used. Pade 
approximants are rational functions of the form 



Rm/n{4>) 



E i=0 °i 
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(11) 



that can be used to approximate functions. These ap- 
proximants typically have better-behaved asymptotics, 
i.e. stay closer to the approximated function, than Tay- 
lor expansions because of their rational structure. An 
extensive expose on Pade approximants can be found in 
Ref. For our study, we will assume 



V{<I>)=R M/N {<t>)., 



(12) 



where again <fi is in units of Mp with (f> presently zero. 
Specifically, we use Pade series Rq/i, Ri/i and i?o/2, as 
these form an exhaustive set of lowest order and next- 
to-lowest order non-trivial expansions with two or three 
parameters. Higher orders are unmotivated given the 
known difficulty for data to constrain more than two dark 
energy /quintessence evolution parameters @, ( as 
will also be evident from our results). 

Pade series have poles, but, as will be discussed in 
the Results Section, data constrains models so that the 
presence of poles is not felt. 

To enable comparison between our results for the two 
different parameterization classes, the priors for the Pade 
series case are set by evaluating the MacLaurin expansion 
of the Pade series, identifying the order coefficients, and 
using the Taylor-series priors for those, i.e. 
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This does not limit us to a finite region, so we additionally 
require |&i| < 2. 
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C. Tracker potentials 

Cosmological tracker potentials/solutions have been 
studied in detail by numerous authors 0, [I?], [H, O H3| • 
These potentials are such that the late-time evolution of 
the field can be essentially independent of initial condi- 
tions, thus providing a possible solution to the coinci- 
dence problem. This behaviour is achieved through a 
type of dynamical attractor solution, and the conditions 
for it to be possible given a particular potential have 
been given and studied in detail by Steinhardt et al. [H| . 
Defining T = V"V/V' 2 , where prime denotes a derivative 
with respect to the field, the two sufficient conditions for 
a potential to possess a tracker solution are 



r > l- 
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(16) 
< 1.(17) 



The first of these conditions ensures convergence to the 
tracker solution (i.e. perturbations away from it are sup- 
pressed), and the second ensures an adiabatic evolution 
of the field that is necessary for the first condition to be 
applicable (and is what one would expect of a function 
that is to maintain a dynamical attractor independent of 
initial conditions). 

If these conditions are fulfilled, the field will eventually 
approach the tracker solution (unless the initial quint- 
essence energy density is too low), and the equation of 
state will then evolve according to 

- 2(r - i) 

W0 ~ ^tracker = ! + 2 (F - 1) ' ^ ' 

possibly breaking away from the tracker solution if ei- 
ther of the conditions later become violated. In assessing 
whether tracking is taking place, one also has to check 
whether the actual evolution on the tracker potential cor- 
responds closely to the tracker solution. An illustration 
of tracker behaviour can be seen in Fig. [T] 

We additionally impose the condition < where 
u>b is the background energy density. This is to ensure 
a possible solution of the coincidence problem by having 
the dark energy density grow with respect to the matter. 
This third condition is usually avoided by specifying the 
tracker condition as T > 1 rather than Eq. ifTB)) . The 
reason for not choosing T > 1 as our condition is related 
to our numerical treatment, and is discussed further in 
Section HVCll 

As we need a non-zero second derivative of the poten- 
tial with respect to the field for T to fulfil the tracker 
conditions, we restrict ourselves to the quadratic poten- 
tial and the Pade series for the tracker viability analysis. 



III. OBSERVABLES 

The observables used are essentially geometric and are 
hence related to the comoving distance for an FRW cos- 
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FIG. 1: Examples of the behaviour of the equation of state 
(here called wq) for a tracker potential. The oscillating curves 
correspond to higher (solid) and slightly lower (dash-dotted) 
initial conditions at high redshift for compared to the 
tracker solution value. The initial velocity at high redshift 
is assumed to be zero. The amplitude of oscillations in w(z) 
around the tracker solution (thin-dashed curve originating at 
wq — 0) decays exponentially with decreasing ln(l + z), and 
the evolution thus approaches the tracker solution regardless 
of the different initial conditions. Although not directly cor- 
responding to our models, the figure illustrates qualitatively 
the tracker property. Reproduced from Ref. [jq|. 



mology described by the parameter vector ©, given by 



where 



E(z; 0) = [n m (l + zf + (1 - n m )e F ^ 



1/2 



and 



F{z;®) = 3 [ (l + i^(z';0))dln(l + z'). 
Jo 



(19) 



(20) 



(21) 



In accordance with our assumptions, these expressions 
assume zero curvature and that quintessence and non- 
relativistic matter are the only relevant components for 
the redshifts we consider. 

We have not included growth-of-structure observa- 
tions, which are not yet competitive with the measures 
we do use (see e.g. Ref. 2l| for a directly-comparable 
example) . 



A. SNIa luminosity— redshift relation 



The luminosity distance is given by 

d L (z;®) = Pl ^ Q ) =(l + z ) r (z;&). (22) 



4 



The apparent magnitude m(z; 0) of a type la supernova 
can be expressed as 



m(z; 0) = M + 5 log 



10 



rf L (s;e) 
Mpc 



25. 



(23) 



where M is the absolute magnitude of SNIa (supposing 
they are standard candles). This can be rewritten as 



m{z; 0) = M + 5 log 10 V L {z; 0) 



(24) 



where M = M - 5 log 10 (if Mpc) + 25 = M — 
51og 10 (/i 70 ) + 43.16 [where h 70 = iJ /(70km/s/Mpc)]. 
Note that some authors define this quantity somewhat 
differently. 

We use the 115 measurements of m(z) mea- 
sured/compiled by the SNLS team [l(|, covering the red- 
shift range z = 0.015 to z = 1.01. The observed magni- 
tudes (indexed by i) are given by 



+ a(si - 1) - (3ci 



(25) 



where m B is the rest-frame B-band magnitude at maxi- 
mum B-band luminosity, and s and c are the shape and 
color parameters. These are derived from the light-curve 
fits and are reported by the SNLS team. The parame- 
ters a and (3 are free parameters and should be varied in 
cosmological fits. However, as they are independent of 
cosmology [22i , we fix them to the SNLS best-fit values 



a = 1.52 + 0.14, 
13 = 1.57 + 0.15, 



(26) 
(27) 



without introducing any bias, and include their uncer- 
tainty in the magnitude uncertainties we use. 

For comparison to our previous paper where the pa- 
rameter r] is used, the parameter M. = M-^ icss — rj, with 
M^ icss the estimate of intrinsic supernova magnitude in 
Riess et al. 



B. CMB peak-shift parameter 



The CMB peak-shift parameter [22 



TZ(z dcc ; 0) 



£l m H a r(z dcc ; 0) 



(28) 



measures an overall linear shift of the CMB power spec- 
trum in multipole space, induced by the effect has 
on the angular-diameter distance to the surface of last 
scattering at z = z dcc . The position of the first power 
spectrum peak is essentially a measure of this distance. 

We use the recent WMAP3 data [1] as analyzed by 
Wang and Mukherjee j24j], who found 



K(z dcc = 1089) = 1.70 ± 0.03 . 



(29) 



C. Baryon acoustic peak 

The standard Big Bang scenario predicts that close to 
the surface of last scattering, baryons and photons act as 
a fluid with acoustic oscillations from the competition be- 
tween gravitational attraction and radiation pressure. As 
the photons decouple, these acoustic oscillations should 
be frozen in the baryon and dark matter distributions. 
One would thus expect an excess of power in the power 
spectrum of luminous matter at a scale corresponding 
to the acoustic scale at last scattering (see e.g. Ref. [25| 
and references therein). Independent first detections of 
this baryon acoustic peak were made by the Sloan Dig- 
ital Sky Survey (SDSS) 9] and the 2dF galaxy redshift 
survey [26] . The SDSS team defined a distance quantity 

A(z BAO ; &) = V^( f° V ^ BAo;0 M , (30) 

which we will use for our analysis. The measurement 
(independent of dark energy model) from the SDSS lu- 
minous red galaxy power spectrum is 

^bao = 0.35) = 0.469 (^) ~°' 35 ± 0.017 , (31) 

which, assuming the WMAP3 mean value n$ = 0.95 [|[, 
yields A{z = 0.35) = 0.474 ± 0.017. 



IV. DATA ANALYSIS 

A. Parameter estimation 

The parameter space we study will be 

= (AL 0o , potential parameters) , (32) 

and we will consistently let D denote the number of free 
parameters in a model. The parameter estimation is car- 
ried out using an MCMC approach, as outlined in our 
previous paper [6j. The posterior probability of the pa- 
rameters 0, given the data and a prior probability dis- 
tribution n(0), is 

P(0|data) = I e -(xsNiJ©)+x^MB(©)+x^o(©))/2 n (0) , 

(33) 

(34) 
(35) 



where 



x|NIa(0) 

Xcmb(©) 
Xbao(O) 



E (mj - m(z l ; 0)) 
^2 



i=l 
(7?-obs 



■R(z dec ;&)Y 



(A obs - A(z BAO ;&)y 



(36) 



Here, we sum over all A^Nia data points for the SNIa 
data, and Z = J £(data|0)n(0)d0 is a normalization 
constant, irrelevant for parameter fitting. Overall, we 
have 115(SNIa) + l(CMB) + l(BAO) data points. 
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B. Model selection 

Separate from the question of parameter estimation is 
the question of parameter necessity, i.e. model selection. 
We again employ an approximate model selection crite- 
rion, the Bayesian Information Criterion (BIC) [27l. |28||. 
given by 



BIC 



-21n£ n 



DlniV, 



(37) 



where £ max is the likelihood of the best-fitting parame- 
ters for that model, D the number of model parameters, 
and N the number of datapoints used in the fit. Models 
are ranked with the lowest value of the BIC indicating 
the preferred model. A difference of two for the BIC 
is regarded as positive evidence, and of six or more as 
strong evidence, against the model with the larger value 
[29L I30I ] . The BIC has also been deployed for dark energy 
model selection in Ref. 1311] . 



C. Tracker viability 

1. Identifying tracker solutions 

To classify general scalar field evolutions as coming 
from a tracker potential capable of solving the coinci- 
dence problem or not, we need to test for both tracker 
conditions and whether the field evolves according to the 
tracker solution. As these conditions are approximate in 
nature, we must specify some e > 0, 5 > such that if 



r > 1 



1 - w h 
6 + 2w h 



(II 



< e . 



din a 

\w<j> — ^tracker] < $, 



(38) 

(39) 

(40) 
(41) 



are all fulfilled for some range in redshift over which we 
require the field to be in the tracker solution, the po- 
tential is classified as a tracker potential. To provide a 
satisfactory solution to the coincidence problem, the field 
should have < % while in the tracker solution. This 
condition is automatically satisfied if the tracker condi- 
tions are fulfilled with T > 1 and the field is in the tracker 
solution. However, in our analysis there is some room for 
fields with > since the field is allowed to deviate 
slightly from the tracker solution, and we also consider 
r > 1 — (1 — w\ ) )/(6 + 2wb) as tracking rather than T > 1 
that is typically used. Cases satisfying the former but not 
the latter are generally disfavoured because they would 
correspond to > in the tracker solution and hence 
not be very successful for solving the coincidence prob- 
lem. In our set-up this is not necessarily true, and this 
is the reason for not choosing the more commonly-used 
latter criterion. Instead, we ensure a solution to the co- 
incidence problem by enforcing < w^. In particular, 



we require T > 5/6 and < since we are concerned 
with the matter-dominated epoch. 

Note that we are not connecting our analysis directly 
with any specific particle physics model and its initial 
conditions at early times, and assessing whether the 
present-time observables are highly insensitive to vari- 
ations in those initial conditions. We are only addressing 
the question whether the (essentially late-time) evolution 
of quintessence is more consistent with such a class of 
tracker potentials, or with a class that does not have such 
behaviour. As the shape of the potential at high redshifts 
is almost unconstrained by data (see also e.g. Ref. [5J]), 
we adopt the viewpoint that a suitable true tracker po- 
tential with insensitivity to initial conditions can always 
be made to coincide with our low-redshift behaviour. 



2. Tracker or Non-Tracker? 

To assess whether models which exhibit tracker solu- 
tion behaviour are favoured by data over models which 
do not, we need some quantity to measure this prefer- 
ence. A well-defined and well-motivated quantity is pro- 
vided within the framework of Bayesian model selection 
[HHlill, where the Bayes factor 



Si 



P(D|Afi) _ P(Mi|D) n(M 2 ) 



P(B\M 2 ) P{M 2 \D) II(Mi) 



(42) 



simply the relative power of Model 1 (Mi ) over Model 2 
(M 2 ) in explaining the observed data D given the prior 
model probabilities IT(Mi) and II(M2), can be used to 
perform this type of comparison. 

For the purposes of assessing the viability of tracker 
solutions for explaining the observed data, we will define 
the following models: 

Mi = {V is a tracker potential} , (43) 
M2 = {V is not a tracker potential} . (44) 

As these two models are disjoint subsets of the model 
space, the Bayes factor can be estimated from Monte 
Carlo Markov chains: letting / pos t be the fraction of 
chain elements from the posterior distribution satisfying 
the tracker criteria, and / pr i the corresponding fraction 
for the prior distribution, the Bayes factor is given by 



B 



12 



/post(l /pri) 
/pri(l /post) 



(45) 



since the fractions of tracker and non-tracker chain el- 
ements must sum to one for both prior and posterior. 
In the limit of equal fractions in prior and posterior, 
B12 = 1, whereas in the limit of complete suppression 
of tracker models in the posterior (so that / pos t = 0) we 
have B12 = in which case Model 2 is infinitely favoured 
over Model 1. 

A standard reference scale for the strength of evidence 
given by the Bayes factor is the Jeffreys scale [29| , shown 
in Table [J 
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ln(Bi 2 ) 


Evidence against Model 2 


0- 1 
1 - 2.5 
2.5-5 

> 5 


Worth only a bare mention 
Positive evidence 
Strong evidence 
Decisive evidence 



TABLE I: The Jeffreys evidence scale. 



We compute the uncertainties in the Bayes factor fol- 
lowing a procedure described in Appendix A. 

The method presented above treats tracker behaviour 
as a Boolean one-parameter property. It is thus insensi- 
tive to intrinsic biases of the combined potential param- 
eterization and parameter priors in fulfilling the different 
tracker criteria, as well as how close to the tracker cri- 
terion limits models typically fall. It would be possible 
to go further and estimate the distributions of param- 
eters measuring each of the three tracker criteria. We 
outline a possible procedure for this in Appendix B, but 
present data do not appear to justify such a sophisticated 
approach and we do not pursue this further here. 



V. RESULTS 
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FIG. 2: One and two-dimensional likelihood distributions for 
a cosmological constant model (A). Solid lines are marginal- 
ized ID likelihoods and dotted lines mean ID likelihoods. 
Solid 2D contours represent 68.3% and 95.4% regions of the 
marginalized distribution, and shading reflects the mean dis- 
tribution. 



A. Parameter estimation 

We present the probability distributions for the fitted 
models in Figs. [2HS1 Marginalized parameter constraints 
and best-fit values are given in Tables ITT1 and IIII1 Plots of 
some dynamical properties of the best-fit models can be 
found in Figs. [7] and [5J The results are discussed further 
below, and model comparison carried out in the following 
subsection. 



1. Cosmological constant (D — 2) 

The probability distributions for the cosmological con- 
stant case are shown in Fig. [2] The parameter constraints 
in Table HTl are improved by roughly a factor of two com- 
pared to our previous analysis [6]. They differ slightly 
from the results of Ref. [33| using the same dataset, al- 
beit within uncertainties. This is most likely due to their 
different treatment of SNLS SNIa errors. 



2. Skater (D = 3) 

The likelihood distributions are shown in Fig. [3] for the 
full dataset, and in Fig. 2]for SNLS alone. Note the sym- 
metry in due to the dependence only on </>q. The 
degeneracy between Vq and 4>o present in our previous 
analysis (where |0o| was positively correlated with Vo) 
is no longer apparent with the full dataset, while still 
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FIG. 3: As Fig. [2] for a 'skater' model, a constant potential 
with kinetic energy. 



being visible if we use supernovae alone. This degen- 
eracy stems from the fact that with supernovae we are 
really only sensitive to an effective quintessence equation 
of state [34:, 35] , which the data require to be close to — 1 . 
Thus, increasing the kinetic energy of the field must be 
compensated by an increase in potential energy to main- 
tain the same effective equation of state. 

Additionally, the mild preference in the Riess et al. 
'gold' data for a non-zero <j>o is not present in the SNLS 
sample, despite the </>o-Vo degeneracy being present. In- 
stead, the likelihood distribution is essentially flat in 
4>o . This could be a reflection of the better qual- 
ity/homogeneity of the SNLS sample over Riess et al. 
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FIG. 4: As Fig.[2]for a constant potential with kinetic energy. 
SNLS data only. Note that here the prior flki n (z > 1) < 0.5 
is applied only up to z — 2. 
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FIG. 5: As Fig. [2] for a linear potential. 



(another possibility is the difference in redshift cover- 
age). In the previous analysis, these two effects conspired 
to give a different best-fit value of Vb in the skater sce- 
nario (Vb = 0.74) compared to the cosmological constant 
(where Vq = = 0.69). That we here do not feel the 
degeneracy is to some degree linked to our prior limit- 
ing f2ki n (z > 1) < 0.5 now being applied to much higher 
redshifts, restricting the range of allowed <fio. However 
the new data do reduce the degeneracy significantly on 
their own (we checked by doing the analysis without the 
prior on f2 k i n ). Also, using only the SNLS data with 
^kin(l < z < 2) < 0.5 (Fig. S]), the flatness of the distri- 
bution in 0o ensures that the best-fit value of Vq in that 
case is only marginally different from that for the full 
analysis, even though the degeneracy is stronger. These 
observations illustrate the need for good-quality data sen- 
sitive to perturbation growth history (e.g. weak lensing) 
to break the 4>q-Vq degeneracy. 



3. Linear potential (D = 4) 

The likelihood distributions are shown in Fig. [51 Note 
the bimodality of the 4>o~Vi distribution, reflecting that 
models are identical under simultaneous change of sign of 
(j)Q and odd-order expansion coefficients. The first change 
from previous constraints [6J] is that the V"o-0o degener- 
acy is now clearly visible in the case of the linear poten- 
tial (there were only hints of it in the previous analysis). 
That is to say, the data quality is getting closer to hitting 
the degeneracy. In addition, we have a degeneracy be- 
tween Vi and 0o , coming from the possibility to achieve 
a particular velocity of the field in the past by either 
changing the present velocity 0o or the slope V\. 

Although not excluding the possibility, the new data 



do not favour a potential where the field is rolling up- 
hill (corresponding to the upper right-hand and lower 
left-hand quadrants of the 0o — V\ distribution). This 
appears to be due to the new SNLS data, which do not 
show a particular preference for a non-zero present field 
velocity, thus not pushing us into these quadrants. It 
would appear that the preference for an uphill rolling 
field found in our previous analysis [|| was an artifact of 
the Riess et al. data. The observational consequences of 
such an uphill rolling field could be interpreted as w < — 1 
if an 'unsuitable' parameterization is used to fit the data 
[34l [36| . It could thus be that the strong w < — 1 prefer- 
ence found in the Riess et al. data (see e.g. Ref. [37|) is 
due to some systematic effect in the data, causing a pref- 
erence for an uphill rolling field and also corresponding 
to a preference for w < — 1 in fits of w(z). This agrees 
with the findings of Nesseris and Perivolaropoulos [371 ] . 
who for three different parameterizations of w find that 
the best-fit w(z) consistently does not cross the phan- 
tom divide line w = — 1 with the SNLS dataset, but does 
with the Riess et al. 'gold' set. The analyses by Barger 
et al. [Hj], Xia et al. [39[ and Jassal et al. [I(| lend sup- 
port to this conclusion as well, as does a recent analysis 
by Nesseris and Perivolaropoulos [41], who however find 
that other cosmological data do gently favour phantom 
divide line crossing provided 0.2 < f2 m < 0.25. 



This also highlights the importance of interpreting 
analyses with care, as we are not probing w(z) directly 
34, 35]. This has been elaborated upon by several au- 
thors in terms of eigenmodes, either as principal compo- 
nents [42] or weight functions [43] ]. 
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FIG. 6: As Fig.[2]for a Pade series with M = 0, N = 1 (7? /i 



4- Pade Ro/i potential (D — 4) 

The likelihood distributions are shown in Fig. [Sj As 
the Ron potential is close to the linear case for small cj>, 

we can use this to compare results. That is, when 0o or 
b\ (which mainly determine the field velocity) are close 
to zero we should expect results to compare well with the 
linear potential which, comparing Fig. [5] with Fig. O we 
see they do. Thus, the discussion above for the linear po- 
tential applies to this case as well. However, as we move 
away from </>o = and b\ — 0, we see that b\ is limited to 
somewhat smaller values than for the linear case (using 
the relation V\ ~ — &061), while the constraints on </>o are 
almost identical. This indicates that data prefer not to 
move very far away from a linear potential. The other 
main feature of the likelihood distributions are bumps 
found in the </>o-fri distributions. These are a feature of 
the likelihood distribution, but the exact size depends on 
our prior enforcing f2ki n (-z > 1) < 0.5 up to high red- 
shifts. 

Pade series, by construction, have poles. One might be 
concerned about how this affects our results if the field 
reaches a pole, but the data is sufficiently constraining 
that the poles are effectively never felt. We tested this 
by doing the analysis with a prior excluding all models 
where a pole is reached before z = 5, and saw no change 
in the results. 



Models with D > 4 



In the three cases with D = 5 (quadratic, Ri/%, and 
^0/2)1 we nn d that the additional parameter is uncon- 
strained by the data and, as in Ref. [6j, we learn nothing 
useful about parameters from these models. Their prin- 
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FIG. 7: Best-fit potentials as a function of the field, and the 
field and potentials as a function of redshift. 



cipal interest lies in model comparison, discussed next, 
where the best-fit found can still be used to assess how 
the models compare in explaining the data. 



B. Model comparison 

The BIC values obtained for all models are shown in 
Tables HI] and IIIII Note that although some parameteri- 
zations have unconstrained parameters, their BIC value 
can be evaluated with Eq. (|37)) from the best fit found in 
our Monte Carlo Markov chains. It is clear that the cos- 
mological constant, showing a BIC difference of at least 
4.6 compared to the other models, is positively favoured. 
This is a strengthening compared to our previous analysis 
where this value was 4.0. In fact, the best-fit x 2 changes 
only marginally between models, thus providing strong 
evidence against linear/Pade Rq/i and higher-order po- 
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TABLE II: Marginalized median and best-fit model parameters and BIC values for the cosmological constant (Qa = Vo/p c ,o) 
and Taylor-series parameterizations. Best-fit values are given in parentheses when differing from the median. Note that the 
likelihood distribution is symmetric under simultaneous change of sign of <f>o and odd-order potential expansion coefficients. 
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TABLE III: Marginalized median and best-fit model parameters and BIC values for the Pade series parameterizations. Best-fit 
values are given in parentheses when differing from the median. 



tentials whose extra parameters add no value. An inter- 
esting feature of the new dataset is that it much more 
strongly disfavours a quadratic potential over the other 
Taylor expansions than just the Riess et al. data. Like- 
wise, the lowest-order Pade expansion is favoured by the 
same amount compared to the higher-order Pade expan- 
sions. 

The best- fit cosmologies (Figs.[7]and[5|) now show more 
convergence in their dynamical properties, although still 
exhibiting increasing variation with redshift. In particu- 
lar, we find that where previously the evolution of 17^ for 
the best-fit quadratic potential was such that O^, stayed 
between 0.75 and 0.96 (for < z < 2), the evolution is 
now very reasonable (see Fig. [5]). The strong evolution 
previously seen in is now more limited, reflecting the 
order-of-magnitude smaller best-fit values for cf>o and V\ 



(though the overall compression of the uncertainties is 
much less than this). 

All best-fit models fall into the 'freezing' category of 
Caldwell and Linder [l^]. For the skater model this 
behaviour is built-in, but it is somewhat intriguing in 
terms of naturalness that the best-fit linear potential ex- 
hibits freezing while at the same time rolling downhill 
(see Figs. [7] and [8|). The potentials with curvature incor- 
porate this best-fit behaviour by making the field reach 
the potential minimum in the recent past (around z = 0.5 
to z = 1), thus providing a braking force to precipitate 
the accelerated expansion of the universe. This situation 
would appear somewhat more natural from a dynamical 
point of view, and it could be that the best-fit linear po- 
tential is trying to approximate this, though data is un- 
able to sufficiently constrain the models with curvature 
in the potential. On the other hand, model selection us- 
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FIG. 8: Best-fit and w$ as a function of redshift. 



ing the BIC also strongly disfavours these models. The 
conclusion must be that complementary or better-quality 
data is needed to resolve this possible contradiction. 

If the linear-potential results stand up, they will 
put the well-motivated models of quintessence based on 
pseudo-Nambu-Goldstone bosons (pNGBs) [44| and simi- 
lar models under pressure, as these rely on a thawing field 
that is becoming dynamical and cosmologically dominant 
in the present epoch. However a field just passing the po- 
tential minimum fits well with the pNGB picture, as well 
as other tracker-type potentials that show a cross-over 
behaviour, such as the SUGRA [45| and Albrecht-Skordis 
[4(j| potentials where the field is starting to feel a curva- 
ture in the potential at late times. These models exhibit 
early quintessence [47|, and can thus be constrained us- 
ing big bang nucleosynthesis and CMB observations [48| . 
It will be interesting to see what future data, including 
those sensitive to perturbation growth and supernovae, 
can tell. 

These observations are in line with studies by 
e.g. Bludman [l9| and Linder [2(|, who both conclude 
that quintessence generically cannot be described by 
slow-roll, and that tracking must break down and move 
towards slow-roll in the recent past (begging the question 
why this is happening precisely now). 



C. Tracker viability 

In carrying out the tracker viability analysis, we con- 
sider four implementations in all by combining two 




FIG. 9: Model average of ln(_Bi2) for tracking required be- 
tween redshift 1 and 10, as a function of e and 8. 



choices of conditions. The first is to demand either that 
the field remains in the tracker regime until the present, 
or that it is allowed to break out of tracking after a red- 
shift of z = 1. The second is to consider two different 
upper limits for the redshift range where the field is re- 
quired to be in the tracker regime, namely z = 2 and 
z = 10; the former more or less represents where the data 
actually lie, while the latter extrapolates the potential to 
higher redshifts. 

We find that all four cases give qualitatively the same 
outcome, and so focus on just one choice, where tracking 
is imposed between z = 10 and z — 1. 

The model average of lni?i2, denoted (I11.B12), for this 
scenario is shown in Fig.[9l for different combinations of e 
and S. For combinations of sufficiently-small e and 5, no 
models satisfying our tracker conditions are found in the 
prior and/or posterior (with those e and 5 limits different 
for the different parameterizations). We exclude these 
cases from our model average, as they effectively corre- 
spond to an infinite uncertainty in the derived value for 
In B12 ■ A very small fraction of the models feel the pres- 
ence of a pole at a redshift lower than the upper tracker 
regime redshift, and are also excluded. We also point out 
that for Pade R0/1, T = 2. Thus, the first two tracker 
conditions are automatically fulfilled, corresponding to 
a delta-function prior on C± and C2 in the language of 
Appendix B. One might consider this a strong bias, and 
hence we exclude this parameterization from our Bayes 
factor model average, and thus use the quadratic, R\n 
and R0/2 potentials to arrive at our conclusions. 

It is clear from Fig. [9] that the average indication 
is in favour of tracker behaviour over non-tracker be- 
haviour. The smallest value of the Bayes factor in the 
figure is 0.98. Limiting our attention to the region where 
e < 0.1,(5 < 0.1, and hence the tracker conditions are 
best obeyed, the smallest value is 2.9. This general trend 
is seen in all four cases we analyze, with the strongest 
preference for tracking in the case presented. However, 
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the model uncertainties in (lni?i2) are comparable to 
(In B12) (particularly for small e and 5) and a firm conclu- 
sion thus cannot be drawn. (As a side note, the Poisson 
uncertainties are relatively small and contribute at most 
on the order of 10% to the total uncertainties.) 

The possible preference for tracker fields is in 
contrast with the commonly-discussed expecta- 
tion > —0.8 for trackers, based on gen- 
eral inverse-power-law series potentials [l8[ (here, 
w c ^ = f a w^(a)fl c f > (a)da/ f a fi^(a)da). While this 
seems to indicate that tracker potentials are disfavoured 
by current data, our results suggest that the data may 
act somewhat more strongly against non-tracker models 
than against tracker ones. 

VI. CONCLUSIONS 

We have updated parameter constraints on the quint- 
essence potential along with cosmological parameters us- 
ing recent SNLS supernova luminosity-redshift data, the 
WMAP3 CMB peak-shift parameter measurement, and 
the SDSS measurement of baryon oscillations. The pre- 
ferred field dynamics appear robust under the different 
parameterizations used. 

We find that, compared to our previous work [6j, pa- 
rameter constraints are improved by roughly a factor of 
two. We also find that linear-potential models where the 
field rolls uphill, although not excluded, no longer pro- 
vide the best fit to the data. The previous mild prefer- 
ence for these models appears to have been an artifact of 
the Riess et al. 'gold' SNIa data. This observation agrees 
with the conclusions of other authors that the SNLS data 
do not particularly favour an equation of state crossing 
the phantom divide line, whereas the Riess et al. data do. 
Although higher-order potentials are not constrained by 
the data, those best-fit potentials exhibit 'cross-over' be- 
haviour, feeling a curvature in the potential in the recent 
past. This qualitatively agrees with some well- motivated 
tracking quintessence models. 

From the point of view of model selection, the cos- 
mological constant is now even more strongly favoured 
compared to the dynamical models we consider (see also 
Refs. [H, H^|). The models with curvature in the po- 
tential are also strongly disfavoured as compared to the 
constant and linear potentials, which appear dynamically 
less natural in the context of the complete evolution ex- 
pected from high redshift. 

We employ a model selection framework to investigate 
whether potentials that exhibit tracker behaviour at in- 
termediate/late times are favoured by data over those 
potentials that do not. We conclude that although our 
results show some indication that tracker behaviour is 
favoured, the model uncertainty on the result is too large 
to draw any firm conclusion. We note that if the dynam- 
ics of our higher-order best-fit potentials and the prefer- 
ence for a tracking potential both stand up in the light 
of new data, the coincidence problem in the context of 



quintessence may simply appear in a new guise — why 
is the field starting to slow-roll now? 

It will be interesting to see how future perturbation 
growth data will help break degeneracies, and, combined 
with supernova and CMB data, constrain quintessence 
models and potentially change the model selection pic- 
ture as well. 
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APPENDIX A: UNCERTAINTY IN TRACKER 
BAYES FACTOR ESTIMATES 

For simplicity of notation we define E = In B12 in this 
Section. The uncertainty in our estimate of E will con- 
sist of two components: Poisson noise from sampling the 
distribution, and model uncertainty. The Poisson noise 
goes as 

CT L = /pri/^pri, (Al) 

CT / post = Jpost/AWt, (A2) 

where JV pr i and N post are the total numbers of samples 
drawn from the prior and posterior distribution respec- 
tively. Accordingly, using standard error propagation 
with Eq. (|45|) . we have that 

4 12 = D 2 a 2 c + C 2 a 2 D + 2B 12 cov(C, D) , (A3) 

with C = /post/ (1 - /post) and D = (1 - / pr i)//pri so that 
B12 = CD. Additionally, we have 

°° - <w£?' (A4) 

°*d = (A5) 

J pri 
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In the absence of knowledge about the covariance be- 
tween C and D, we can place an upper limit on the Pois- 
son uncertainty, 

>2 



v 2 Bl2 < {Da c + Ca D y 



(A6) 



We use this upper limit as our estimate for the Poisson 
uncertainty. The corresponding uncertainty in E is then 



a E = 



C D 



+ 



/post(l /post) /pri(l /pri 



(A7) 



The model average of E over M models is given by 
(note that this quantity is denoted (hiB-12) in the main 
text) 



E = 



M 



with an associated uncertainty 



M(M - 1) 



(A8) 



(A9) 



We will now have an 'error on the error' from the Poisson 
uncertainty, given by 



EM 



y M{M - 1) J2i(E — Ei) 2 
so our final estimate of E will be 
E = E 



(A10) 



(All) 



± 



! E,(E-E t y 

M(M — 1) 



M(M — 1) J2i(E — Ei 



APPENDIX B: TRACKER PROBABILITY 
DISTRIBUTIONS 

Here we briefly describe a possible extension of the 
tracker analysis carried out in this paper, though we be- 
lieve application to present data would be premature. 

To address the model uncertainty in the Bayes fac- 
tor model average, we consider the probability distribu- 
tions of the parameters that determine whether a model 
is classed as a tracker. In more detail, we can define three 
different 'tracker functions' 



Ci[z tr ] 

C 2 [z tr ] 
C 3 [z tr ] 



min (r(z) — 5/6) 

zGz tr 



T(z) 



-idT(z) 



dim 



(Bl) 
(B2) 



max 

z6z tr 



max \w<j,{z) ~ «WkerO?)| , (B3) 

ZGZ tr 



where z tr is the redshift range for which the field is re- 
quired to exhibit tracker behaviour, and record their val- 
ues for all elements in our MCMC chains. Note that 



we do not include a function corresponding to the con- 
straint ua < 0, as maxw^(z) will be a function of C\ and 
C3. From this we obtain the posterior probability distri- 
bution P(Ci, C 2 , C3|II(©)) given the prior distribution 
11(0) for our primary cosmological parameters 0. Run- 
ning the MCMC for the prior distribution as well, we 
obtain the prior distribution n(Ci, C 2 , C3|n(0)). 

We are then in a position to do importance sampling 
(see e.g. Appendix B in Ref. [5(| for a brief introduction) 
using the prior and posterior we have calculated. We 
can change priors for C\ , C2, C3 from those induced by 
n(0) to whichever we like and obtain the corresponding 
new posterior distribution, since we only need to divide 
out the prior distribution and multiply by the prior of our 
choice (with the exception of parts of parameter space cut 
out by the primary prior n(0) or very poorly sampled). 
A potential problem with this approach is that optimal 
sampling of the posterior distribution in Ci, C2, C3 is not 
necessarily achieved by optimal sampling in the primary 
parameters, and sufficient statistics may take a long time, 
i.e. many chain elements, to accumulate. 

Setting natural priors for these new parameters may 
be perceived as difficult (although not manifestly more 
arbitrary than for other phenomcnological parameteriza- 
tions). A simple way of setting the priors is to argue 
that we should be equally likely to draw a parameter 
value that fulfils the corresponding tracker criterion, as 
one that doesn't. For instance, if we assume Gaussian 
priors, we get 



P{C 2 ) 
P(C 3 ) 



1 



27T(TCi 

2 



2 



2tt<t c , 



exp 



■ exp 



• exp 



C? 



2^ 
CI 



2^3 



9{C 2 



0(C 3 



(B4) 



(B5) 



(B6) 



where 6 is the Heaviside step function (C 2 and C3 are re- 
stricted to non- negative values by definition). The stan- 
dard deviations <rc 2 an d f c 3 are set by then demanding 

/ P{C 2 )AC 2 - / P(C 2 )dC 2 , (B7) 

JC 2 <e Jc 2 >£ 

[ P(C 3 )dC 3 = f P(C 3 )dC 3 . (B8) 

JC,<& JCa>6 



The case of C± is different, since we only have one in- 
equality to fulfil (r > 5/6). Hence, we need to put a 
cut-off at some value to determine the standard devia- 
tion. One could of course assign, for example, flat priors 
in the same fashion. 

Using this method, we can thus obtain a posterior 
distribution P(C\ 1 C 2 ,C 3 ) for a given prior distribution 
n(Ci, C2, C3) of our choice, thus allowing a removal 
of correlation biases intrinsic to particular parameteri- 
zations, which should reduce model uncertainty. This 
method allows us to perform parameter estimation on 
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Ci, (72,6*3 as well as model selection by calculating the 
Bayesian evidence. It is of course applicable to gen- 
eral dynamical cosmological properties one might wish 
to study. Carrying this out in practice can however be 
involved since we might not be sampling efficiently in 



the MCMC, and performing model selection in a robust 
manner would require specialized code to address the 
sampling inefficiency and to handle the use of a binned 
distribution. 
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